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1 Introduction 

Ferromagnetism occurs when Mn is randomly substituted for more than 
about 2 percent of the cations of several III-V compound semiconductors. 
Although only a few host materials have been explored at present, this prop- 
erty is likely shared by most III-V semiconductors. In this Chapter we will 
discuss some of the theoretical pictures that are being developed to explain 
the magnetic and transport properties of these materials. Our development 
will be based on a phenomenological model that has been used with great 
success to explain the sensitivity of bulk and layered (II,Mn)VI semiconduc- 
tor optical properties to external magnetic fields. (Ferromagnetism does not 
occur for Mn in undoped II- VI hosts.) The low energy degrees of freedom 
in this model are holes in the semiconductor valence band and one S = 5/2 
local moment for each Mn ion. 

Interest in these ferromagnets was heightened by the demonstration sev- 
eral years ago that ferromagnetic transition temperatures [0 in excess of 
100 K can be achieved in (Ga,Mn)As. It has been further heightened recently 
both by the demonstration of long spin-coherence times in semiconductors 

and by the dramatic and rapid development of new information storage 
technology based on magnetotransport effects in ferromagnetic metals |^ . It 
seems clear that semiconductors have many potential advantages over met- 
als for devices based on the magnetotransport effects that occur in itinerant 
electron ferromagnets, principally because they present a wider canvas for 
their creative manipulation by some combination of impurities, gates and 
optical excitation. It is likely that important applications for these materials 
will be found only if ferromagnetism at room temperature can be achieved. 
The very recent discovery [Q of ferromagnetism at temperatures close to 
1000 K in (Ga,Mn)N has fueled hopes that these materials will indeed have 
technological impact. Our focus here, however, is on the physics of these fer- 
romagnets; we therefore concentrate mainly on the properties of (Ga,Mn)As 
and (In,Mn)As which have been studied most extensively ||^. 

It is generally accepted that Mn acts as an acceptor when it substitutes 
for a cation in a III-V semiconductor lattice, leaving a Mn^+ ion which has a 
half-filled d-shell with angular momentum L = and spin S = 5/2. It is also 
generally accepted that ferromagnetism occurs in these materials because of 
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interactions between Mn local moments that are mediated by holes in the 
semiconductor valence band. There is, therefore, a lot of similarity between 
the ferromagnetism of these materials and that of lanthanides and actinides 
and their compounds in which f-electron moments are coupled by d-band 
itinerant electrons. There are also similarities between these materials and 
the manganite compounds that have been extensively studied |^ in recent 
years in part because of the large increase in resistance that occurs when T 
exceeds the Curie temperature, the so-called colossal magnetoresistance ef- 
fect. In (III,Mn)V ferromagnets, however, the local moments appear on only a 
small fraction of the atomic sites arranged randomly. In addition the itinerant 
electron density is also low, even lower than the Mn local moment density. As 
we explain later this property is likely important in selecting ferromagnetic 
over glassy magnetic order. The participation of itinerant electrons in the fer- 
romagnetism of these diluted magnetic semiconductors (DMS) adds to their 
richness, leading in particular to electronic transport properties that are very 
sensitive to the magnetic state of the material. The physics of ferromagnetic 
semiconductors is in a sense intermediate between that of rare earth magnets 
and that of manganites in that the spin-splitting of itinerant electron bands 
due to their exchange coupling with local moments is comparable to their 
Fermi energies, rather than being much smaller than band Fermi energies as 
in the rare earth case or larger as in the manganite case. 

The Chapter is organized as follows. In Section ^ we briefly summarize 
the main experimental properties of (III,Mn)V ferromagnets. Several different 
but related approaches that have been explored in an effort to gain insight into 
these materials are outlined in Section ^. The remaining sections of this re- 
view Chapter deal with the development of the semi-phenomenological model 
we favor in which the low energy degrees of freedom are exchange-coupled 
valence-band holes and S = 5/2 Mn local moments that carry a negative 
charge. The simplest version of this model is one in which the randomly 
distributed Mn ions are replaced by a uniform continuum, thus completely 
neglecting disorder. In Section ^ we discuss physical predictions based on 
a mean-field treatment of this disorder-free model and demonstrate that it 
successfully describes a number of non-trivial properties of (Ga,Mn)As and 
(In,Mn)As ferromagnets, including their anomalous Hall conductivities. In 
Section |^ we discuss collective excitations of these ferromagnets within the 
disorder-free model, demonstrating that the simple mean-field-theory is rea- 
sonably reliable for typical parameters of current samples but must fail at 
large carrier densities and also in the limit of very strong exchange coupling. 
In Section || we discuss the results of Monte Carlo calculations that describe 
the effect of collective fluctuations of Mn moment orientations. The method 
can deal with some of the complications and additional physics, including the 
possibility of non-coUinear ground states, that enters when disorder is added 
to the theoretical model. A brief summary is given in Section 
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In Sections ^, j^, and || we have been able to present only a small fraction 
of theoretical data obtained using our approach. An extensive survey of pre- 
dicted physical properties of DMS's will be available at http : // unixl2 . fzu.cz/ 
ms/index.php web-pages, launched by the authors of this Chapter in collab- 
oration with Jan Kucera, Byounghak Lee, and Jairo Sinova. 

2 Properties of (III,Mn)V Ferromagnets 

In this section we discuss some important properties of (III,Mn)V ferromag- 
nets that have been established by current experiments. Thorough recent 
reviews of the properties of these materials have been prepared by Dietl, 
Matsukura, and Ohno @J^. Our objective in this section is to summarize 
the observations that are most important in constraining theoretical descrip- 
tions. There is at present considerable activity related to the growth and 
characterization of these materials. We expect rapid progress to be achieved 
in the near future in exploring the range of possible behaviors and relating 
them more precisely to molecular-beam-epitaxy growth and post-growth an- 
nealing protocols. We list below a number of properties that appear to be 
safely established. A theoretical picture that is able to explain most of these 
properties is outlined in the following sections. 

• Electron paramagnetic resonance and optical experiments P,p^ demon- 
strate that S = 5/2 local moments occur for dilute concentrations of 
Mn in GaAs. These experiments demonstrate that the Mn local moment 
model is correct. 

• The Mn-induced states near the Fermi energy play a key role in the origin 
of ferromagnetism and in the magnetotransport properties of (III,Mn)V 
DMS's. According to photoemission studies pT| , p"2|jT^ , those states have 
As 4p character, i.e., can be associated with the host semiconductor 
valence band states. The angle-resolved photoemission experiment also 
showed a negligible shift in the heavy- and light-hole bands with Mn 
concentration x < 7%. 

• Ferromagnetism is not observed for Mn concentrations smaller than ~ 
0.01 1^. This property demonstrates that ferromagnetism does not oc- 
cur when all valence-band holes are trapped on individual Mn ions or on 
other defects. Antisite defects, for example, are common in semiconduc- 
tor samples grown by low-temperature molecular beam epitaxy (MBE). 
For very dilute Mn concentrations, electron spin resonance experiments 
demonstrate that most holes arc trapped not at the Mn acceptors, but 
at other defects. 

• The ground-state magnetization, M{T — 0), per Mn ion can exceed 

for larger values of a: ||^ . Since the magnetization contribution from anti- 
ferromagnetically coupled valence band holes tends to partially compen- 
sate the Mn local moment magnetization, ferromagnets with these large 
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values of M{T = 0) likely have ground states with (nearly) fully aligned 
Mn local moments. 

• To date the largest ferromagnetic transition temperatures in (Ga,Mn)As 
occur for x ~ 5%. The current record is Tc ^ llOK The drop in 
critical temperatures at higher x values may be related to Mn clustering 
or may have a more fundamental origin. There is a correlation between 
large values of M(T = 0) and high T^'s. 

• (III,Mn)V thin film ferromagnets grown under compressive strain have 
their magnetic easy axis in the plane, while ferromagnets grown under 
tensile strain have their magnetic easy axis in the growth direction 
External magnetic fields ^ 100 mT are sufficient to align the magneti- 
zation along the hard axis 0,0. These properties can be explained by 
well understood strain effects in the spin-orbit coupled valence bands and 
demonstrate that the macroscopic properties of these ferromagnets are 
sensitive to details of the valence band electronic structure. 

• The ferromagnetic critical temperature and the temperature dependence 
of the magnetization are altered by post-growth annealing and are sensi- 
tive to the details of the annealing protocol Q . 

• These ferromagnets have large anomalous Hall resistivities [|], demon- 
strating that the itinerant valence bands are full participants in the mag- 
netism. The large anomalous Hall resistivities reflect the strong spin-orbit 
coupling that is present at the top of the valence band in zincblende semi- 
conductors. 

• The semiconductor valence bands are spin-split in the ferromagnetic 
state. These semiconductor ferromagnets exhibit strong magnetoresis- 
tance effects, like tunnel magnetoresistance that are characteristic 
of itinerant electron ferromagnets, again demonstrating that the itinerant 
electrons are full participants in the magnetism. 



3 Theoretical Approaches 

Ferromagnetism is a collective effect due to interactions between electrons. If 
it were possible to do so, we would explain its microscopic origins in a par- 
ticular class of materials by solving the many-electron Schrodinger equation 
directly, given the position of all the nuclei. Fortunately this uninteresting 
direct approach is impossible, now and forever more, because of the macro- 
scopic number of interacting electronic degrees of freedom. Instead we must 
resort to some combination of approximation and phenomenological mod- 
eling, settling on the correct approach only after careful comparison with 
experiment. This is the art of condensed matter science; an intricate tango 
between theory and experiment leading to a conclusion that cannot be antic- 
ipated while the dance is in progress. In some cases, fractional quantum Hall 
systems and possibly cuprate superconductors for example, the low-energy 
degrees of freedom in terms of which observable physical properties are best 
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described are not simply related to the bare electrons which appear in the 
many-electron Hamiltonian. Often, however, the low-energy degrees of free- 
dom are more obvious. 

A practical approach to many-electron physics that is often successful is 
spin-density-functional (SDF) theory, in which many-body effects appear in 
exchange-correlation potential contributions to effective independent-particle 
Hamiltonians. SDF theory has the advantage that it is a first principles ap- 
proach without any phenomenological parameters. Among its many achieve- 
ments is a generally satisfactory description of itinerant electron ferromag- 
netism in transition metals. SDF theory has been applied to (III,Mn)V 
ferromagnetism by Sanvito et al. and by van Schilfgaarde and Mryasov. To 
date these calculations have been performed using the local density approx- 
imation (LDA) of SDF theory, an approximation that is not reliable when 
local moments are formed, i.e., when strong correlations suppress fluctua- 
tions in the number of electrons in the d-shell or the f-shell of a particular 
atom. LDA-SDF theory supercell jl^ calculations predict that majority spin 
d-electrons of a Mn atom substituted on a cation site of GaAs lie at the 
Fermi energy, rather than lying well below the Fermi energy as they would 
if the half- filled d-shell formed a 5 = 5/2 local moment. Similar results have 
been obtained in coherent potential approximation band-structure calcula- 
tions JTrt for (In,Mn)As, and are clearly in disagreement with experiment 
in both cases. In local moment systems, it is necessary to account for the 
increase in instantaneous cite energy when the occupancy of a localized or- 
bital is increased. The LDA-I-U method has been developed to mitigate this 
deficiency of SDF theory and has recentlylQ been applied to (III,Mn)V fer- 
romagnets and finds a dominant As 4p orbital weight at the Fermi energy, 
consistent with the photoemisson experiment. 

It appears likely that a lot of detailed information on the electronic prop- 
erties of (III,Mn)V ferromagnetic semiconductors can be reliably obtained 
from LDA-I-U SDF calculations and we expect this approach will play an 
important role in the modeling of these materials in the future. Our work, 
however, follows a semi-phenomenological strategy, starting from a model 
in which the local-moment character of the Mn d-orbitals is asserted rather 
than derived. These models are adapted from ones used to describe 

the optical properties of (II,Mn)VI semiconductors. The low-energy degrees 
of freedom in the kinetic-exchange model ||2l| , ^ we employ are 5 = 5/2 local 
moments representing half-filled Mn d-shells, and holes in the Mn valence 
band (see Fig. |^). Since Mn^+ ions should act as acceptors when substituted 
on the cation sites of III-V semiconductors, we would expect one hole per 
Mn if no other charged defects were present in the system. However, antisite 
defects are common when III-V semiconductor films are grown by MBE at 
the low temperatures required to prevent Mn segregation. The hole density 
and the Mn density are therefore taken as separate sample-dependent quan- 
tities, to be determined experimentally. The Hamiltonian of this model is 
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specified in detail in the following paragraph. There has also been theoretical 
work on these materials based on a still simpler model where holes are 
assumed to hop only between Mn acceptor sites, where they interact with 
the Mn moments via phenomenological exchange interactions. These models 
have some advantages in getting at the physics of the dilute Mn limit, and 
can also easily be adapted to include the holes that are localized on ionized 
antisite defects rather than Mn acceptors 




Fig. 1. Model for (III,Mn)V semiconductors: local magnetic moments (Mn'^"'") with 
spin 5 = 5/2 are antiferromagnetically coupled to itinerant carriers (holes) with 
spin s = 1/2. 



Like any phenomenological model, the one we use is defined most fun- 
damentally by its low-energy degrees of freedom. Also important, however, 
is the Hamiltonian that acts in the implied Hilbert space. The length scales 
associated with holes in these compounds are still long enough that a k ■ p, 
envelope function, description |^ of the semiconductor valence bands is ap- 
propriate and we take that approach here. The operators in terms of which 
the phenomenological Hamiltonian is expressed include the spin operator Sj 
for the 5 = 5/2 local moment on site / and the multi-band envelope function 
hole spin density operator s(r). The following key terms are included in the 
minimal version of the model Hamiltonian: 

a) The coupling of the Mn spin to the external magnetic field, gfiB X)/ '^i ' 

b) The band Hamiltonian of the host III-V semiconductor, usually described 
using a multi-band envelope function formalism p5| , |26[ |. For many prop- 
erties it is necessary to incorporate spin-orbit coupling in a realistic way. 
Six- or eight-band models that include the 'split-off' band and/or the 
conduction band are sometimes desirable. Unlike the local moment mod- 
els, the Hilbert space includes all host lattice sites for each hole rather 
than only localized orbitals centered on the Mn sites. This band Hamil- 
tonian should include the strain effects due to lattice matching between 
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the epitaxially grown (III,Mn)V films and the substrate on which they 
are grown. 

c) Antiferromagnetic exchange coupHng between the Mn^"'' spin and valence- 
band holes, Jpd^iSj ■ s{Ri). This interaction represents virtual cou- 
pling to states that have been integrated out of the model's Hilbert space, 
ones in which electrons are exchanged between the Mn ion d shells and 
the valence band |^,^. The exchange interactions are isotropic to a 
good approximation because the Mn^"*" ion has total angular momentum 
L = 0. Experimental estimates for Jpd vary from 150 ± 40 meV nm^ 

to 68 ± 10 meV nm'^ to 55 ± 10 meV nm'^ [0. Recent experimental 
work limits Jpd to a value toward the lower end of this range and fixes 
its value within perhaps 20%. 

The terms d)-f), listed below, are necessary to describe the crossover to 
the localized limit in which holes are bound either to Mn acceptors or to other 
defects. Note that simpler, impurity-band models assume that the system is 
in this limit from the outset, much as our model assumes from the outset 
that the Mn d-shells form local moments. There are sometimes technical 
difficulties in describing the localized limit with our higher-level model, so 
that considerable simplification arises from using the impurity-band model. 
There is however, a penalty to pay since the model does not apply to the 
regimes of greatest interest in which the band electrons are not localized. 
Even when impurity models do apply, it is difficult to guess at appropriate 
distribution functions for the inter-site hopping parameters that play a key 
role. 

d) The attractive Coulomb interaction between the ionized Mn^+ acceptor 
and a valence-band hole. In an envelope function formalism, central-cell 
corrections to the interaction are necessary to capture the isolated bound- 
acceptor limit accurately 

e) The repulsive Coulomb interaction among holes. This interaction is key 
in screening the ionized Mn^^ acceptors and cannot be neglected ex- 
cept in the completely-localized-hole limit. When it is included, it usu- 
ally must be approximated in a way which avoids artificial hole-hole self- 
interactions. 

f ) The repulsive interaction between holes and ionized antisite (group-V el- 
ement on group-Ill site) defects. The antisite defects compensate for the 
Mn acceptors and reduce the overall hole density, in addition to providing 
an important additional scattering center. In the dilute Mn limit, experi- 
ment suggests that most Mn^+ ions do not have bound holes, possibly 
due to this compensation. When all holes are strongly localized, most Mn 
local moments will be free and the system will not have ferromagnetic 
order. 



The following terms are potentially important in some circumstances. 
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g) The scalar scattering potential that represents the energy difference be- 
tween a valence-band p electron on a host site and a valence-band p 
electron on a Mn site. This effect has normally been excluded since its 
size and sign is not yet know. 

h) Direct exchange interactions between Mn ions on neighboring sites. These 
terms result from microscopic processes in which exchange of electrons 
between the valence band and two nearby Mn d-shells is correlated. Terms 
of this type are known to be important in (II,Mn)VI semiconductors, but 
appear to be less important in (III,Mn)V semiconductors. 

i) Direct coupling of band electrons to external magnetic field. 

In the rest of this Chapter we discuss only pictures of (III,Mn)V fcrro- 
magnetism that follow from the minimal model that includes only the a)-c) 
Hamiltonian terms. 

4 Mean-field-theory predictions 

In this and the following section we make an important approximation that 
achieves a drastic simplification. We will refer to this as the continuum Mn 
approximation, although it has sometimes been referred to as a virtual crystal 
approximation. It is motivated by the observation that the Fermi wavelength 
of the valence-band electrons is typically longer than the distance between 
Mn ions, mainly because the Mn acceptors are compensated and also partially 
because there are four occupied valence bands. When the Mn ion distribution 
is replaced by a continuum with the same spin-density, randomness is com- 
pletely eliminated from the minimal model. This approximation has many 
elements in common with the dynamic mean-field-theory (DMFT) approxi- 
mation, applied pO| to these ferromagnets recently by Chattopadhyay et ai, 
although the DMFT does retain some of the consequences of randomness 
neglected here and in the following section. The possibility of starting with 
a description based on an approximation where the random Mn distribution 
is replaced by a continuum emphasizes an essential difference between ferro- 
magnetic semiconductors and classical spin glass systems in which dilute Mn 
local moments are distributed randomly in a metallic host. In both cases it 
is true that the interaction between local moments mediated by the itiner- 
ant electrons is oscillatory and ferromagnetic for separations smaller than the 
itinerant electron Fermi wavelength. However, the Fermi wavelength is longer 
than the distance between local moments in the doped semiconductor case, 
whereas it is shorter in the spin glass case . Each Mn ion interacts ferro- 
magnetically with several of its neighbors. The continuum Mn approximation 
will fail in the limit of dilute Mn ions and also if the exchange interaction 
between band and Mn spins is too strong. This approximation does not al- 
low for the sensitivity of magnetic properties to annealing protocols that has 
been established in experiment. It does, however, seem to be reliable in the 
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limit of principle interest, that of high Mn densities and high critical temper- 
atures, where the holes are metallic and their interaction with Mn acceptors 
will be effectively screened. The merit of this approximation is that it enables 
quantitative prediction of many physical properties. In this section we dis- 
cuss three important properties, the ferromagnetic transition temperature, 
the magnetic anisotropy energy, and the anomalous Hall conductance. This 
work described below is motivated by the view in science, it is ultimately 
up to experiment to decide on the reliability of any approximation made in 
modeling a physical system. As we will point out, the utility of the contin- 
uum Mn approximation is strongly supported by observations. In fact, it is 
not a surprise that this approximation is a good starting point, given its 
success in describing the influence of external fields on the properties of the 
closely related paramagnetic (II,Mn)VI semiconductors [ pO|p^ . The present 
section makes in addition a mean-field approximation by ignoring correlations 
between Mn and band spin configurations. 

Our mean-field theory is derived in the spin-density-functional frame- 
work and leads to a set of physically transparent coupled equations . The 
effective magnetic field seen by localized magnetic ions consists of an exter- 
nal magnetic field and the mean kinetic-exchange-coupling contribution from 
spin-polarized carriers, 

HesiRl) = ffext + Jpd{s{Rl))/9HB , (1) 

where {s{Rj)) is the carrier spin density at Mn sites, and g is the g-factor of 
the local moments. The mean spin polarization of a magnetic ion is given by 

{S)i = ~SBs{SgfiBH,s{Ri)/kBT)H,fi{Ri) , (2) 

where Bs{x) is the Brillouin function and HcsiRi) is the unit vector along 
the direction of the effective magnetic field defined in Eq. (Q). The itinerant- 
hole spin density is determined by solving the Schrodinger equation for holes 
which experience a kinetic-exchange effective Zeeman field h{r). The field 
h(r) is non-zero only in the ferromagnetic state and, in the continuum limit, 
reads 

h{r) = NMn{r) {S){r) , (3) 

where Nmu — 4a;/a^^ is the Mn density in Mna;IIIi_^V zincblende semicon- 
ductors with a lattice constant aic- For inhomogeneous systems such as quan- 
tum wells or superlattices, the itinerant holes experience also an electrostatic 
potential due to heterostructure confinement and external bias (if present). 
Using the local-spin-density approximation (LSDA), itinerant hole-hole in- 
teraction can be accounted for by including an additional spin-dependent 
one-particle potential in the Schrodinger equation. 
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4.1 Ferromagnetic transition temperature 

In homogeneous DMS systems, the hole-spin density (s) and the kinetic- 
exchange potential h are related at small h by 

(''^)=-T^^h. (4) 

Here, g* is the hole g-factor and Xf is the interacting hole magnetic suscep- 
tibility, 



where Etot/V is the total energy density of the itinerant-hole system. The 
Curie-Weiss transition temperature, obtained from Eqs. (|l]) - (||) in the 
Hext = limit, is 

ksTc = — r2 . (6) 

To understand the qualitative physics implicit in this Tc-equation (^, 
we discuss first the magnetic susceptibility expressions of a model itinerant 
electron system with a single spin-split band and an effective mass m*. The 
kinetic-energy contribution E^"^ to the total energy gives 

where fci? is the Fermi wavevector. The exchange energy of the spin-polarized 
parabolic-band model adds a contribution 

dh^ ATT^eh^ ' 

where e is the dielectric constant of the host semiconductor. At high hole 
densities p, the kinetic-energy term dominates and is proportional to the 
Fermi wavevector, i.e., to p^^^. Equations ^ and (||) also show that the 
band contribution to the mean-field Tc increases linearly with m* while the 
exchange enhancement of Tc is proportional to (m*)^. Note that correlation 
effects, not discussed here in detail, suppress the mean-field Tc by only ~ 1% 
for typical experimental hole densities {p ~ 0.1 nm~^) in bulk (III,Mn)As 
ferromagnets. 

To obtain quantitative predictions for Tc, it is necessary to evaluate the 
kinetic and exchange contribution to the itinerant hole susceptibility using a 
realistic six-band Kohn-Luttinger model instead of the parabolic band 
model. The Hamiltonian contains the spin-orbit splitting parameter Ago and 
three other phenomenological parameters, 71, 72, and 73, whose values for 
the specific III-V host can be found, e.g., in Refs. Results (7|j3|,||,|3 
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Fig. 2. The band (kinetic energy) contribution to the mean-field ferromagnetic 
critical temperature Tc for Mn concentration x = 5% is plotted as a function of 
hole density p for InAs, GaAs, AlAs, and GaN host semiconductors. 

are plotted in Fig. || as a function of the hole density for InAs, GaAs, AlAs, 
and GaN host semiconductors. In the density range considered, only the two 
heavy-hole and two light-hole bands are occupied in the arsenides. However, 
the mixing between these four bands and the two spin-orbit split-off bands 
is strong and must be accounted for. In GaN, spin-orbit coupling is weak 
and all six bands are occupied by holes. The numerical data are consistent 
with the qualitative analysis based on the parabolic band model: the band 
(kinetic-energy) contribution to Tc follows roughly the p^^^ dependence, the 
exchange enhancement of ~ 10% is only weakly density dependent. The Tc 
values at a given density are ordered according to the heavy-hole and light- 
hole masses in the arsenide hosts. For GaN, with all six bands occupied, the 
simple model of a parabolic spin-split band is less instructive. Yet the large 
numerical Tc's in this material are consistent with the large heavy-hole mass, 
nearly twice as large as in AlAs. 
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The mean-field prediction for the critical temperature agrees quantita- 
tively with the experimental value of 110 K measured in Mn-doped GaAs 
with Mn concentration x = 5% and p — 0.35 nm~'^. Thermal fluctuations ne- 
glected by the mean-field theory, discussed in the following section, reduce the 
theoretical Tc estimate by less than 5% ||3^ , explaining the quantitative suc- 
cess of the mean-field theory in this sample. The same analysis finds approx- 
imately a Tc suppression of approximately 20% compared to mean-field 
theory due for (Ga,Mn)N, implying that room temperature ferromagnetism 
may occur in III-V DMS. 



4.2 Magnetic anisotropy 

Experiments in (III,Mn)V DMS's have demonstrated that these ferro- 

magnets have remarkably square hysteresis loops and that the magnetic easy 
axis is dependent on epitaxial growth lattice-matching strains. The physi- 
cal origin [^J2^ of the anisotropy energy in our model is spin-orbit coupling 
in the valence band. Even in mean-field theory, we find that the magnetic 
anisotropy physics of these materials is rich and that easy axis reorienta- 
tions can occur as a function of sample parameters including hole density or 
epitaxial growth lattice-matching strains. 

Magnetic anisotropy in the absence of strain is well described by a cubic 
harmonic expansion truncated at sixth order, an approximation commonly 
used in the literature | ]39| on magnetic materials. The corresponding cubic 
harmonic expansion for total energy of a system of non-interacting holes in 
the presence of the effective field h is 



= ^totiiioo)) + + hihi + hihi) + KT him 



where h is the unit vector along the field h. The cubic anisotropy coefficients 
K"^ and ifj^ are related to total energies for h along the high symmetry 
crystal directions by following expressions: 

.^ca_ 4[i^tot((110))-i?tot((100))] 

^ V 
^ca ^ 27£:tot((lll)) - 36i;tot((110)) + 9igtot((100)) ^ ^^^^ 

MBE growth techniques produce (III,Mn)V films whose lattices are locked 
to those of their substrates. X-ray diffraction studies Q have established that 
the resulting strains are not relaxed by dislocations or other defects, even for 
thick films. Strains in the (III,Mn)V film break the cubic symmetry assumed 
in Eq. (^). However, the influence of MBE growth lattice-matching strains 
on the hole bands of cubic semiconductors is well understood [p5| and we 
can use the same formal mean-field theory as in the previous subsection to 
account for strain effects on magnetic anisotropy. 
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Fig. 3. Cubic magnetic anisotropy coefficients Kl'^ and K2'^ as a function of hole 
density p. 

We turn now to a series of illustrative calculations intended to closely 
model the ground state of (Ga,Mn)As. For Mn density A^Mn = 1 nni~^ 
{x « 5%), h K, 140 meV at zero temperature. This value of h is not so 
much smaller than the spin-orbit splitting parameter in GaAs [^,0 (Ziso = 
341 meV), so that accurate calculations require the six-band Luttinger model 

. Even with A^Mn fixed, our calculations show that the magnetic anisotropy 
of (III,Mn)V ferromagnets is strongly dependent on both hole density and 
strain. The hole density can be varied by changing growth conditions or by 
adding other dopants to the material, and strain in a (Ga,Mn)As film can be 
altered by changing substrates. The cubic anisotropy coefficients (in units of 
energy per volume) for strain-free material are plotted as a function of hole 
density in Fig. |^. The easy axis is nearly always determined by the leading 
cubic anisotropy coefficient ATJ^*^, except near values oip where this coefficient 
vanishes. As a consequence, the easy axis in strain free samples is almost al- 
ways either along one of the cube edge directions (ATJ'* > 0), or along one 
of the cube diagonal directions [Kf^ < 0). Transitions in which the easy 
axis moves between these two directions occur twice over the range of hole 
densities studied. (Similar transitions occur as a function of h, and therefore 
temperature, for fixed hole density.) Near the hole density p = 0.01 nm~'^, 
both anisotropy coefficients nearly vanish and a fine-tuned nearly perfect 
isotropy is achieved. The slopes of the anisotropy coefficient curves vary as 
the number of occupied bands increases from 1 to 4 with increasing hole den- 
sity. This behavior is clearly seen from the correlation between oscillations of 
the anisotropy coefficients and onsets of higher band occupations. 

Six-band model Fermi surfaces are illustrated in Figs. ^ and ^ by plotting 
their intersections with the kz — plane at p = 0.1 nm^'^ for the cases of 
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(100) and (110) ordered moment orientations. The dependence of quasipar- 
ticle band structure on ordered moment orientation, apparent in comparing 
these figures, should lead to large anisotropic magnetoresistance effects in 
(Ga,Mn)As ferromagnets. We also note that in the case of cube edge orienta- 
tions, the Fermi surfaces of different bands intersect. This property could have 
important implications for the decay of long-wavelength collective modes. 



0.10 




-0.10 I ' • ' • ' • ' 1 

-0.10 -0.05 0.00 0.05 0.10 

Fig. 4. Six-band model Fermi surface intersections with the fcz = plane for p — 
0.1 nm~'^ and h = 140 meV. This figure is for magnetization orientation along the 
(100) direction. 



In Fig. 1^ we present mean- field theory predictions for the strain-dependence 
of the anisotropy energy at /i = 140 meV and hole density p — 0.35 nm~'^. 
According to our calculations, the easy axes in the absence of strain are along 
the cube edges in this case. The relevant value of the in-plane strain produced 
by the substrate-film lattice mismatch, 

eo = -, 11 

a/ 

depends on the substrate on which the epitaxial (Ga,Mn)As film is grown. 
The most important conclusion from Fig. ^ is that strains as small as 1% 
are sufficient to completely alter the magnetic anisotropy energy landscape. 
For example for (Ga,Mn)As on GaAs, eo = -0.0028 at x = 0.05. The 
anisotropy has a relatively strong uniaxial contribution, even for this rel- 
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Fig. 5. Six-band model Fermi surface intersections with the kz = plane for the 
parameters of Fig. ^ and magnetization orientation along the (110) direction. 




-0.10 -0.05 0.00 0.05 0.10 
eo 

Fig. 6. Energy differences among (001), (100), (HO), and (111) magnetization ori- 
entations vs. in-plane strain eo at h = 140 meV and p = 0.35 nm^'^. For compressive 
strains (eo < 0), the system has an easy magnetic plane perpendicular to the growth 
direction. For tensile strains (eo > 0), the anisotropy is easy-axis with the preferred 
magnetization orientation along the growth direction. The anisotropy changes sign 
at large tensile strain. 
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atively modest compressive strain, which favors in-plane moment orienta- 
tions 0:^^, in agreement with experiment ||^. A relatively small (^1 kJ 
m~'^) residual in-plane anisotropy remains which favors (110) over (100). 
For X — 0.05 (Ga,Mn)As on a x = 0.15 (In,Ga)As buffer the strain is ten- 
sile, eo = 0.0077, and we predict a substantial uniaxial contribution to the 
anisotropy energy which favors growth direction orientations , again in 
agreement with experiment For the tensile case, the anisotropy energy 
changes more dramatically than for compressive strains due to the depopu- 
lation of higher subbands. At large tensile strains, the sign of the anisotropy 
changes, emphasizing the subtlety of these effects and the latitude which 
exists for strain-engineering of magnetic properties. 



4.3 Anomalous Hall effect 

The mean-field description of hole bands in the presence of exchange cou- 
pling to the localized Mn moments provides a starting point for building a 
theory of transport in (III,Mn)V ferromagnets. Here we concentrate on the 
anomalous Hall effect which is an important sample characterization tool 
in ferromagnetic systems. The Hall resistivity of ferromagnets has an ordi- 
nary contribution, proportional to the external magnetic-field strength, and 
an anomalous contribution usually assumed to be proportional to the sam- 



ple magnetization. In our approach |40 , the anomalous Hall conductance of 



a homogeneous ferromagnet is related to the Berry phase of the electronic 
wavefunction acquired by a cyclic evolution along the Fermi surface. 

In the standard model of the AHE in metals, skew-scattering ^l[] and 
side-jump scattering give rise to contributions to the Hall resistivity 
proportional to the diagonal resistivity p and respectively, with the latter 
process tending to dominate in alloys because p is larger. Our evaluation of 



the AHE in (III,Mn)V ferromagnets is based on a theory |43 of semiclassical 
wave-packet dynamics which implies a contribution to the Hall conductivity 
that is independent of the kinetic-equation scattering term. The interest in 
this contribution is motivated in part by practical considerations, since our 
current understanding of (III,Mn)V ferromagnets is not sufficient to permit 
confident modeling of quasiparticle scattering. The relation of our approach 
to standard theory is reminiscent of disagreements between Smit | |4l[ | and 
Luttinger Q] that occurred early in the development of AHE theory and 
do not appear to have ever been fully resolved. We follow Luttinger in 
taking the view that there is a contribution to the AHE due to the change in 
wavepacket group velocity that occurs when an electric field is applied to a 
ferromagnet. The electron group velocity correction is conveniently evaluated 
using expressions derived by Sundaram and Niu : 

de 

- ^ + {e/h)E X fl. (12) 

The first term on the right-hand-side of Eq. ( p^ is the standard Bloch band 
group velocity. Our anomalous Hall conductivity is due to the second term. 
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proportional to the Berry curvature 17, defined below. It follows from sym- 
metry considerations that for a cubic semiconductor under lattice-matching 
strains and with m aligned by external fields along the (001) growth direction, 
only the z-component of fl is nonzero: 

l?.(n,fc)=2Im[(|^||^)]. (13) 

Here |m„) is the periodic part of the n-th Bloch band wavefunction with the 
mean-field spin-splitting term included in the Hamiltonian. The anomalous 
Hall conductivity that results from this velocity correction is 



e 



2 



dk 



o-AH = j T^fn,kn4n,k) , (14) 

where fn,k is the equilibrium Fermi occupation factor for the band quasi- 
particles. We have taken the convention that a positive ctah means that the 
anomalous Hall current is in the same direction as the normal Hall current. 

This Berry phase contribution to the anomalous Hall conductance occurs 
in any itinerant electron ferromagnet with spin-orbit coupling. To assess its 
importance for (HI,Mn)V ferromagnets we first explore a simplified model 
that yields parabolic dispersions for the two heavy-hole and two light-hole 
bands, and neglect coupling to the split-off band by assuming a large spin- 
orbit coupling p5| , p6t . Detailed numerical simulations accounting for the mix- 
ing of the spin-orbit split-off bands and warping of the occupied heavy-hole 
and light-hole bands ||2|,||| in the (In,Mn)As and (Ga,Mn)As samples 
will follow. Within the 4-band spherical model, the spin operator s — j/3, 
and the Hamiltonian for holes in III-V host semiconductors can be written 
as 



r 5 

Ho = Tr- ill + x72)fc' - 272(fc • jf 
2m 2 



(15) 



where j is the total angular momentum operator and 71 and 72 are the Lut- 
tinger parameters | p5| , ^ . In the unpolarized case (h — 0), the total Hamilto- 
nian, H = Hq — hjz/'i (the external magnetic field is assumed to be in the +z 
direction), is diagonalized by spinors where, e.g., jj,=j-k = ±3/2 for the 
two degenerate heavy-hole bands with the effective mass mhh = "1/(71 — 272). 
The corresponding Berry phase, J (fikn{zt3/2, k) — ±3/2(cos6'fe — 1), is 
largest at the equator (cos0fc = kz/k^i^ — 0) and vanishes at the poles 
(|cos0fc| = 1) of the spherical Fermi surface of radius /chh. Because of the 
band degeneracy, the anomalous Hall conductivity (|lj) vanishes in the h = 
limit. The effective Zeeman coupling present in the ferromagnetic state both 
modifies the Fermi surface shapes and renormalizes the Berry phases. Up to 
linear order in h we obtain that k^^ = fchh ± hm\-ih/ (2h^kYih) cos6k and that 
the Berry phase is reduced (enhanced) by a factor [1 ^ 2mft,/(972?i^A:j^jj)]. A 
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similar analysis for the light-hole bands leads to the total net contribution to 
the AHE from the four bands whose lower and upper bounds are: 



47r2?i' 



-h{3Tr^p) ^/^TOhh < CTAH < 



-/l(37r2p)-l/322/3^j^j^ 



(16) 



Here p = k^^^/^Tr^ (1 + y/ "^ih/^hh) is the total hole density and mih = 
m/(7i + 272) is the light-hole effective mass. The lower bound in Eq. ( p^ 
is obtained assuming mih ^ Whii while the upper bound is reached when 
mih « TOhh- 

Based on the above analysis we draw the following conclusions: The 
anomalous velocity due to the Berry phase can have a sizable effect on the 
AHE in (ni,Mn)V ferromagnets. The solid line in Fig. shows our ana- 
lytic results for the GaAs effective masses mhh = 0.5me and mih — O.OSme. 
Note that in experiment, anomalous Hall conductances are in order of 1-10 
fl~^ cm~^ and the effective exchange field h ^ 10— 100 meV. A large ctah is 
expected in systems with large heavy-hole effective mass and with the ratio 
wih/whh close to unity. 
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Fig. 7. Illustrative calculations of the anomalous Hall conductance as a func- 
tion of the band-splitting effective Zeeman field for hole density p — 0.35 nm^^. 
The dotted-dashed curve was obtained assuming infinitely large spin-orbit coupling 
and the decrease of theoretical ctah with decreasing spin-orbit coupling strength is 
demonstrated for A^o ~ 1 eV (dashed line) and A^o = 341 meV (solid line). 



So far we have discussed the limits of infinitely strong spin-orbit coupling 
and weak effective exchange field, relative to the hole Fermi energy. In the 
opposite limits of zero spin-orbit coupling or large h, ctah vanishes. This 
implies that the anomalous Hall conductivity is generally nonlinear in the 
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exchange field or the magnetization. To explore the intermediate regime we 
diagonalized the six-band Luttinger Hamiltonian numerically ^5 26 with the 
spin-orbit gap Ago = 1 eV as well as for the GaAs value Ago = 341 meV. 
Results shown in Fig. |7| confirm that a smaller ctah is expected in systems 
with smaller Ago and suggest that both positive and negative signs of ctah 
can occur, in general. The curves in Fig. ^ are obtained by neglecting band 
warping in III-V semiconductor compounds. The fact that valence bands in 
these materials are typically strongly non-parabolic, even in the absence of 
the field h and in the large Ago limit, is accurately captured by introducing 
the third phenomenological Luttinger parameter 73 |^,^. Numerical data 
including all Luttinger parameters indicates that warping tends to lead to an 
increase of a ah, as seen when comparing solid curves in Fig. |^ and in the top 
panel of Fig. |[ The hole-density dependence of ctah, illustrated in Fig. ||, 
is qualitatively consistent with the spherical model prediction ( |l6[). Also in 
accord with the outlined chemical trends, numerical data in Fig. ^ suggest 
large positive AHE in (Al,Mn)As, intermediate positive ctah in (Ga,Mn)As, 
and a relatively weaker AHE in (In,Mn)As with the sign of ctah that may 
depend on the detail structure of the sample. 

We make now a comparison between our ctah calculations and experi- 
mental data in the (In,Mn)As and (Ga,Mn)As samples, analyzed in detail by 
Ohno and coworkers psPJst . The nominal Mn densities in the two measured 
systems are A^Mn — 0.23 nm^'^ for the InAs host and A^Mn — 1-1 nm^'^ for the 
GaAs host, yielding saturation values of the effective field h — 25±3 meV 
and h = 122 ± 14 meV, respectively. The low-temperature hole density of the 
(Ga,Mn)As sample, p = 0.35 nm~'^, was unambiguously determined ||] from 
the ordinary Hall coefficient measured at high magnetic fields of 22-27 T. 
Since a similar experiment has not been reported for the (In,Mn)As sample 
we estimated the hole density, p — 0.1 nm~^, by matching the measured 
ferromagnetic transition temperature Tc = 7.5 K to the density dependent 
mean-field Tc- 

As demonstrated in Fig. ^, our theory explains the order of magnitude 
difference between experimental AHE in the two samples (ctah ~ 1 cm~^ 
in (In,Mn)As and ctah ~ 14 il^^ cm~^ in (Ga,Mn)As). The calculations are 
also consistent with the observed positive sign and monotonic dependence of 
Ctah on sample magnetizations [D . 

We take the agreement in both magnitude and sign of the AHE as a strong 
indication that the anomalous velocity contribution dominates the AHE in 
homogeneous (HI,Mn)V ferromagnets. This Berry phase term, which is in- 
dependent of quasiparticle scatterers, is relatively easily evaluated with high 
accuracy, enhancing the utility of the Hall measurement in sample charac- 
terization. The success of this model also supports the use of the simple 
mean-field approximation discussed in this section, in which Mn ions are rep- 
resented by a uniform density continuum, to describe at least the ground 
state of these ferromagnets. 
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Fig. 8. Full numerical simulations of ctah for GaAs host (top panel), InAs host 
(bottom panel), and AlAs host (inset) with hole densities p = 0.1 nm~^ (dotted 
lines), p — 0.2 nm"'^ (dashed lines), and p = 0.35 nm~ (solid lines). The filled 
circles in the top and bottom panels represent measured AHE values. The 

saturation mean-field h values for the two points were estimated from nominal 
sample parameters |^,^ . Horizontal error bars correspond to the experimental un- 
certainty of the Jpd coupling constant. The measured hole density in the (Ga,Mn)As 
sample is p = 0.35 nm~^; for (In,Mn)As, p = 0.1 nm~^ was determined indirectly 
from the sample's transition temperature. 



5 Collective excitations within a continuum picture 

5.1 Beyond mean-field theory and RKKY interaction 

The power and the success of the mean-field picture employed in the previous 
section lies in the fact that it is a simple theoretical approach which makes 
it easy to calculate many observable quantities numerically. Mean-field the- 
ory, however, neglects correlation between local-moment spin configurations 
and the free-carrier state and, therefore, fails to describe the existence of 
low-energy long-wavelength spin excitations, among other things. Because of 
its neglect of collective magnetization fluctuations, mean-field theory, e.g., 
always overestimates the ferromagnetic critical temperature. There are many 
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examples in itinerant electron systems where mean-field theory overestimates 
ferromagnetic transition temperatures by more than an order of magnitude 
and it is not a priori obvious that mean-field theory will be successful in 
(III,Mn)V ferromagnets. Indeed, we will find that the multi-band character 
of the semiconductor valence band plays an essential role in enabling high 
ferromagnetic transition temperatures in these materials. 

In this section we identify the elementary spin excitations, determine their 
dispersion, and discuss implications for the Curie temperature 1 46 47| , ^ , ^9|j5C[| . 
The starting point of our analysis is the itinerant-carrier-mediated ferromag- 
netic interaction between local magnetic moments. Such an interaction is 
provided by the familiar Ruderman-Kittel-Kasuya-Yoshida (RKKY) theory. 
The RKKY picture, however, only applies as long as the perturbation induced 
by the Mn spins on the itinerant carriers is small. As we will derive below, the 
proper condition is A <^ ep where A = NunJpdS is the (zero-temperature) 
spin-splitting gap of the itinerant carriers due to an average effective field in- 
duced by the Mn ions, and is the Fermi energy. This condition is, however, 
never satisfied in (III,Mn)V ferromagnets, partially because the valence-band 
carrier concentration p is usually much smaller than the Mn impurity density 
A^Mn- A related drawback of the RKKY picture is that it assumes an instan- 
taneous static interaction between the magnetic ions, i.e., the dynamics of 
the free carriers are neglected. We will see below that this dynamics is im- 
portant to obtain all types of elementary spin excitations. As a consequence, 
RKKY theory does not provide a proper description of the ordered state in 
ferromagnetic DMSs. 

As in the previous section, we use here the minimal model including terms 
a)-c) of Section ^ and employ the Mn continuum approximation. Extensions 
to the minimal model may be important in some circumstances. They are, 
however, not essential for the general discussion in the present section, which 
will attempt to explain the considerations that determine when collective 
effects neglected by mean-field-theory are important. 



5.2 Independent spin-v^rave theory for parabolic bands 

The main idea of our theory is to derive an effective description for the Mn 
spin system by integrating out the valence-band carriers and to look for fluc- 
tuations of the Mn spins around their spontaneous mean-fleld magnetization 
direction (which we choose as the z-axis). Using the Holstein-Primakoff (HP) 
representation | |5l| , we express the Mn spins in terms of bosonic degrees of 
freedom. We expand the effective action up to quadratic order, i.e., we treat 
the spin excitations as noninteracting Bose particles. From the corresponding 
propagator we deduce the dispersion of all elementary spin excitations. 

To keep the discussion transparent we start with a two-band model for the 
itinerant carriers with quadratic dispersion. Later, in Section 5.5, we extend 
our theory to a model with a more realistic band structure described by a 
six-band Kohn-Luttinger Hamiltonian. 
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For small fluctuations around the mean-field magnetization, we can write 
the spin operators as 



S+{r) « b{r)y'2NMnS (17) 
S-{r)^b\r)^2NMnS (18) 
S'{r)^NMnS^bHr)b{r) (19) 

with bosonic fields 6^(r), b{r). The state with fully polarized Mn spins (along 
the z-direction) corresponds, in the HP boson language, to the vacuum with 
no bosons. The creation of a HP boson reduces the magnetic quantum number 
by one. 

The partition function Z can be expressed as a coherent-state path in- 
tegral in imaginary time over the HP bosons and the valence-band carri- 
ers, which are fermions. Since the Hamiltonian is bilinear in the fermionic 
fields, we can integrate out the itinerant carriers and arrive at an effective 
description in terms of the impurity spin degree of freedom labeled by the 
complex number coherent state labels for the boson fields, z and z. We get 
Z — J Dlzz] exp(— S'off [zz]) with the effective action 

S'cffN = Sbp[zz] - Indet [{G^'^)-^ + SQ-Hzz)] , (20) 

where S'bp[zz] = jj^ dr J d^r zdrZ is the usual Berry's phase term. In Eq. (po|), 
we have already split the total kernel G^^ into a mean-field part {G^^)~^ 
and a fluctuating part SG~^, 

(G^'")--! - {dr - Ai) % + {i\Ho\j) + NunJpdSs!^ (21) 
<5G-.i(zz) = ^ [(zs-. + ZS+) V^NunS ~ 2zzsl^ (22) 

where ^ denotes the chemical potential, and i and j range over a complete set 
of hole-band states (i.e., here, for the model with two parabolic bands, i and j 
label band wavevectors and spin, t;i)j ^nd and sf^ are matrix elements of 
the itinerant-carrier spin matrices. The combination A = Nyi^J^^S defines 
the mean-field energy to flip the spin of an itinerant carrier. The physics of 
the itinerant carriers is embedded in the effective action of the magnetic ions. 
It is responsible for the retarded and non-local character of the interactions 
between magnetic ions. 

So far we have made no approximations. The independent spin- wave the- 
ory is obtained by expanding Eq. ( ^o| ) up to quadratic order in z and z, i.e., 
spin excitations are treated as noninteracting HP bosons. This is a good ap- 
proximation at low temperatures, where the number of spin excitations per 
Mn site is small. 

We obtain (in the imaginary time Matsubara and coordinate Fourier rep- 
resentation) an action that is the sum of the temperature-dependent mean- 
field contribution and a fluctuation action. The latter is 

5eff [zz] = ^ ^ z{k,Vrn)D^^ik,iy„i)z{k,iy„i). (23) 

\k\<.k]j ,m 
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A Debye cutoff ko witli fc|, = 67r^A^Mn ensures tliat we include the correct 
number of magnetic-ion degrees of freedom, |fc| < ku. The kernel of the 
quadratic action defines the inverse of the spin-wave propagator, 



where ^ = (p| — p^)/p is the fractional free-carrier spin polarization, and 
e-f_|(<7) is the energy of spin-up and spin-down valence-band holes, e|,|(q) = 
e,'± Z\/2, and = h^q"^ /{2m*). The second term of Eq. (H) is the the 
energy for a Mn spin excitation in mean- field-theory, J7^^ = = xA. 

It differs from the itinerant-carrier spin splitting by the ratio of the spin 
densities x = p^/(2A^MnS'), which is always much smaller than 1 in (III,Mn)V 
ferromagnets. Mean-field theory is, thus, recovered by dropping the last term 
in Eq. (p^). It is this term that describes the response of the free-carrier 
system to changes in the magnetic-ion configuration. 



5.3 Elementary spin excitations 

We obtain the spectral density of the spin-fluctuation propagator by ana- 
lytical continuation, iiy'm ^ f2 + jO+ and A{k,f2) — Im D{k, n)/TT. In the 
following we consider the case of zero temperature, T = 0. We find three 
different types of spin excitations . 



Goldstone-mode spin waves. Our model has a gapless Goldstone-mode 
branch reflecting the spontaneous breaking of spin-rotational symmetry. The 
dispersion of this low-energy mode for four different valence-band carrier 
concentrations p is shown in Fig. ^ (solid lines). At large momenta, fc — > oo, 
the spin-wave energy approaches the mean- field resuh i?^^^ ^ [2^'^ (short- 
dashed lines in Fig.|9|). Expansion of the T — propagator for small momenta 
yields for the collective modes dispersion. 



where is the Fermi energy of the majority-spin band. In strong and weak- 
coupling limits, Z\ ^ ef and Z\ <C ef, respectively, Eq. (^ ) simplifies to 



fi['^^- -ek + 0(k^) for /l>eF, (26) 



,(1) _ 

.(1) P , 



We note that the dependence of the spin-wave energy on the system param- 
eters, namely the exchange interaction strength Jpd, hole concentration p, 
local-impurity density iVMn, and effective mass m* is different in these two 
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Fig. 9. Spin- wave dispersion (solid lines) for Jpd = Q.06eVnm , m* = Q.5me, 
A^'mh = lnm~^, and four different itinerant-carrier concentrations p — 0.01 nm~"^, 
0.035 nm"^, 0.1 nm"^, and 0.35 nm"^. The ratio A/ep is 2.79, 1.21, 0.67, and 0.35, 
which yields the fractional free-carrier spin polarization ^ as 1, 1, 0.69, and 0.31. 
The short wavelength limit is the mean-field result f2^^ = xA (short-dashed lines) , 
and the long-dashed lines are the result obtained from an RKKY picture. 



limits, indicating that the microscopic character of the gapless collective exci- 
tations differs qualitatively in the two limits. The energy of long-wavelength 
spin waves is determined by a competition between exchange and kinetic 
energies. To understand this in more detail one can impose the spin configu- 
ration of a static spin wave on the Mn spin system, evaluate the ground-state 
energy of the itinerant-carrier system in the presence of the generated ex- 
change field, and compare this with the ground-state energy of a uniformly 
polarized state. The results of this calculation are explained briefly below; 
for details see Ref. [Q. Given the Mn spin configuration, the valence-band 
carriers can either follow the spatial dependence of the Mn spin density in 
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order to minimize the exchange energy, as they do in the strong-couphng 
hmit A ^ ep, or minimize the kinetic energy by forming a state with a ho- 
mogeneous spin polarization , as they do in the weak-couphng hmit A <ti ep- 
The corresponding energy scales are provided by A and ep, i.e., the crossover 
from one regime to the other is governed by the ratio A/ep. 

Stoner continuum. We observe that the frequency Vm is not only present 
in the first term of Eq. (|2^), it enters the third term, too. This is the reason 
why, in addition to the Goldstone mode, other spin excitations can appear 
in our model. They are absent in a static-limit description, i.e., when the 
frequency dependence of the third term of Eq. (p^ ) is neglected. 
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Fig. 10. Stoner excitations and optical spin-wave mode in the free-carrier system 
for Jpd ~ 0.06eVnm^, m* = 0.5me, Nmu = lnm~'^, and p — Q.35nm~^. In an 
RKKY picture these modes are absent. 

We find a continuum of Stoner spin-flip particle- hole excitations. They 
correspond to flipping a single spin in the itinerant-carrier system and, since 
x 1, occur in this simple model at much larger energies near the itinerant- 
carrier spin-splitting gap A (see Fig. |lo| ). For A > ep and zero temperature, 
all these excitations carry spin = +1, i.e., increase the spin polarization. 
They therefore turn up at negative frequencies in the boson propagator we 
study. When A < ep, excitations with both = +1 and = —1 contribute 
to the spectral function. The continuum lies between the curves — Zi — ± 
'2^/ek£F and for A < ep also between —A -t- ± 2y^ ek{ep — /\). 



26 J. Konig, J. Schliemann, T. Jungwirth, and A.H. MacDonald 



Optical spin waves. Wc find additional collective modes analogous to the 
optical spin waves in a ferrimagnet. Their dispersion lies below the Stoner 
continuum (see Fig. [l^) . At small momenta the dispersion is 

The finite spectral weight at negative frequencies indicates that, because of 
quantum fluctuations, the ground state is not fully spin polarized. 

5.4 Comparison to RKKY and to the mean-field picture 

For comparison we evaluate the T = magnon dispersion assuming an RKKY 
interaction between magnetic ions. This approximation results from our the- 
ory if we neglect spin polarization in the itinerant carriers and evaluate the 
static limit of the resulting spin- wave propagator defined in Eq. (|2j). The 
Stoner excitations and optical spin waves shown in Fig. ^ are then not 
present and the Goldstone-mode dispersion is incorrect except when A <ti ep, 
as depicted in Fig. || (long-dashed lines). 

In the mean-field picture, correlations among the Mn spins are neglected. 
The mean-field theory can be obtained in our approach by taking the Ising 
limit, i.e., replacing S ■ s hy S^s^ . As mentioned before, this amounts to 
dropping the last term in Eq. ( p^ ) . The energy of an impurity-spin excitation 
is then dispersionless, = xA (short-dashed line in Fig. ||), and always 

larger than the real spin- wave energy. 

5.5 Spin-wave dispersion for realistic bands 

For a quantitative analysis jsot of the spin-wave dispersion we extend our 
parabolic-band model to a six-band Kohn Luttinger Hamiltonian. The effec- 
tive action for the HP bosons describing the Mn impurity spins is given by 
the same formal expression Eq. (|2^) with the contributions Eqs. (|l]) and ( p2| ) 
to the kernel. The difference is that for each Bloch wavevector i and j now 
label the states in a six-dimensional Hilbert space (instead of two dimensions 
for spin up and down), iJp is the Kohn-Luttinger Hamiltonian, and sf^ and 
s,^- are 6x6 matrices. 

The next step is again an expansion of the effective action up to quadratic 
order in z and z. In the two-band model, where spin is a good quantum 
number, only z(k, iym)z{k, the combinations appear, see Eq. (p3|). Since 
the coherent state labels can can be viewed as boson creation and annihilation 
operators, these contributions are diagonal in total boson number. In the 
presence of spin-orbit coupling, however, spin is no longer a good quantum 
number, and the combinations z{k, Vm)z{—k, —Vm) and z{k, Vm)z{—k, —Vm) 
which increase or decrease the number of HP bosons, come into play. 

Since our aim here is to derive the dispersion relations of the low-energy 
spin waves, rather than to address the full excitation spectrum including the 
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Stoner continuum and the optical spin waves, we take the static limit as 
discussed in the context of the two-band model. After a Bogoliubov transfor- 
mation, we obtain for the spin-wave energy 



iTfc _ 
A 

with the definition 




El 



q Q/3 



/[6^(g)]~/[eff(g + fc)] ^ 



(29) 



(30) 



for cr, cr' = ±. The indices a and /3 label the single-particle eigenstates for 
valence-band carriers at a given wavevector q and q + k, and — (a|s^|/3). 
The remaining task is to evaluate the fractional itinerant-carrier polarization 
^ and the quantities Et~ and Et^ numerically. 
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Fig. 11. Main panel: Spin-wave dispersion for the 6-band model for itinerant-carrier 
density p — 0.35 nm~^, impurity-spin concentration Nmu ~ 1.0 nm~'^ and exchange 
coupling Jpd = 0.068 eV nm~^. Inset: Spin- wave dispersion on a log-log plot (circles) 
and the parabolic fit (solid line). 



In Fig. ^ we show the spin-wave dispersion for wavevectors k along the 
easy axis obtained using parameters valid for (Ga,Mn)As ^j. We observe 
that the effect of E^^ in Eq. (|2^) is negligibly small and can, therefore, 
be dropped. Furthermore, we find that the dispersion is fairly independent 
of its wavevector direction, a property that is usually implicitly assumed in 
micromagnetic descriptions of magnetic materials. 
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Spin stiffness. The quantized energy of a long-wavelength spin wave in a 
ferromagnet with uniaxial anisotropy can be written as 

where K is the anisotropy energy constant, and A denotes the spin stiffness 
or exchange constant. While the anisotropy constant can be obtained from 
the mean-field energy for different magnetization orientations (see previous 
section) , the virtue of the spin- wave calculation is to extract the spin stiffness 
as well. 

In Fig. |T2|we show the spin stiffness A as a function of the itinerant-carrier 
density for two values of Jpj for both the isotropic two-band and the full six- 
band model. We find that the spin stiffness is much larger for the six-band 
calculation than for the two-band model. Furthermore, for the chosen range 
of itinerant-carrier densities the trend is different: in the two-band model the 
exchange constant decreases with increasing density, while for the six-band 
description we observe an increase with a subsequent saturation. 
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Fig. 12. Exchange constant yl as a function of itinerant-carrier density p for the 
six-band and the two-band model for two different values of Jpd = 0.068 eVnm~^ 
(solid lines) and 0.136 eV nm~"^ (dashed lines). The impurity-spin concentration is 
chosen as A^Mn ~ 1.0 nm~^, which yields A — 0.17 eV (solid lines) and A = 0.34 eV 
(dashed lines), respectively. 



To understand this behavior we recall that the two-band model predicts 
a different dependence of A on p in the strong and weak-coupling limits with 
a crossover near A ^ ep, see Eqs. ( p6[ ) and ( p7|). The difference in the trends 
seen for the two- and six-band model in Fig. is explained in part by the 
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observation that, at given itinerant-carrier concentration p, the Fermi energy 
ej? is much smaher when the six-band model is employed, where more bands 
are available for the carriers, than in the two-band case. Furthermore, we 
emphasize that, even in the limit of low carrier concentration, it is not only 
the (heavy-hole) mass of the lowest band which is important for the spin 
stiffness. Instead, a collective state in which the spins of the itinerant carriers 
follow the spatial variation of a Mn spin- wave configuration will involve the 
light-hole band, too. Our calculations show that accounting for the presence 
of this second more dispersive band is essential to understanding the suc- 
cess of mean-field theory. Crudely, the large mass heavy hole band dominates 
the spin-susceptibility and enables local magnetic order at high temperatures, 
while the dispersive light hole band dominates the spin stiffness and enables 
long range magnetic order. The multi-band character of the semiconductor 
valence plays an essential role in the ferromagnetism of these materials. 

5.6 Limits on the Curie temperature 

Isotropic ferromagnets have spin- wave Goldstone collective modes whose en- 
ergies vanish at long wavelengths. 



where k is the wavevector of the mode. Spin-orbit coupling breaks rotational 
symmetry which leads to a finite gap, see Eq. According to our numer- 
ical studies, though, this gap is negligibly small as far as the suppression of 
ferromagnetism by collective spin excitations is concerned and can, therefore, 
be dropped for the present discussion. Each spin-wave excitation reduces the 
total spin of the ferromagnetic state by 1. The coefficient D — 2A/ {NMnS) 
is proportional to the spin stiffness A. These collective excitations are not 
accounted for in the mean-field approximation. If the spin stiffness is small, 
they will dominate the suppression of the magnetization at all finite tem- 
peratures and limit the critical temperature. In this case, the typical local 
valence-band carrier polarization remains finite above the critical tempera- 
ture. Ferromagnetism disappears only because of the loss of long-range spatial 
coherence. 

A rough upper bound on the critical temperature T™'' can be obtained by 
the following argument which accounts for the role of collective fluctuations 
The magnetization vanishes at the temperature where the number of 
excited spin waves equals the total spin of the ground state. 



= Dk^ + Oik^) , 



(32) 




(33) 



where n{f2k) is the Bose occupation number and the Debye cutoff, k^ = 
(GTr^iVivin)"'^^^, ensures the correct number of magnetic ion degrees of free- 
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dom. We therefore find that the critical temperature of a ferromagnet cannot 
exceed 

fcsT™" = ^^Dkl (34) 

for 5 > 5/2 where D is the T = spin-stiffness. To obtain this equation, we 
have assumed that the spin waves can be approximated as non-interacting 
Bose particles, replaced the dispersion by the long- wavelength limit Eq. (p2[), 
and noted that the critical temperature estimate is proportional to Dkjj, 
justifying the use of the classical expression for the mode occupation number 
n/c « ksT / fik — 1/2. These considerations set an upper bound on the critical 
temperature which is proportional to the spin stiffness, a bound not respected 
by mean-field theory. 

To get a qualitative but transparent picture we employ the two-band 
model with parabolic bands, and deduce the spin stiffness from Eqs. ( p6| ) 
and ( p7|) for the strong and weak-coupling regime, respectively. For strong 
coupling, A/ep 3> 1, the exchange coupling completely polarizes the valence- 
band electrons, and we find (using p <^ iN-^i^S) the bound 

rpCoU.S 25* + 1 / p \^ , . 

For small Zi/e_F, the weak-coupling or RKKY regime, exchange coupling is a 
weak perturbation on the band system. In this regime we get 

rricoll.RKKY _ rpMF 25* + 1 / A^Mii ^ ^ 

12(5 + 1)^1 P ) ' ^-^^^ 

i.e., mean-field theory is reliable only for p/NMn ^ 1, as expected since in 
this case the RKKY interaction has a range which is long compared to the 
distance between Mn spins. 

We expect that the qualitative picture derived from the two-band model 
will persist for the six-band model. The actual values of the boundaries be- 
tween the regimes indicated in Fig. ^ will, however, be shifted because of 
the important differences in the microscopic physics that determines the spin 
stiffness of the two models discussed above. We expect that the Tc esti- 
mates derived in the preceding paragraph will be directly applicable to n- 
type carrier-mediated ferromagnets. As a consequence we expect that it will 
be impossible to achieve large ferromagnetic transition temperatures in n-type 
semiconductors with carrier mediated ferromagnetism. The ~ 5% reduction of 
(Ga,Mn)As mean- field Tc due to spin fluctuations, mentioned in Section 4.1, 
was obtained using the six-band model and solving self-consistently Eq. ( ^4| ) 
and I?(r™'i) = D{T = 0)(5)(r™")/'S'- Larger reductions compared to mean- 
field-theory estimates are expected in some hosts, but mean-field-theory re- 
tains a qualitative validity. 
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Fig. 13. Critical-temperature-limit regimes for the two-band model. In the mean- 
field regime Tc is limited by individual Mn spin fluctuations. In the collective 
regimes, the critical temperature is limited by long-wavelength fluctuations with 
a stiffness proportional to the bandwidth for weak (RKKY) exchange coupling and 
inversely proportional to the bandwidth for strong exchange coupling. At the solid 
line T^^ = T^". Dashed lines: expansions for large and small A/ep, Eqs. ( p^ ) and 
(ji^), and the crossover from the RKKY to the strong coupling collective regime. 

6 Collective fluctuations beyond spin wave theory and 
continuum approximation 

In the preceding sections we described the magnetic properties of Mn-doped 
semiconductors by a model of itinerant carriers which are exchange-coupled to 
localized magnetic moments formed by the dopants. An important property, 
resulting from the growth process of such materials, is that the Mn acceptors 
are distributed at random on the cation sites of the underlying crystal lat- 
tice. To this point in this chapter, we have used a continuum approximation 
for the Mn ion distribution that yields a disorder-free problem. The contin- 
uum approximation makes it possible to obtain some results analytically, and 
also crucially simplifies numerical calculations. However, it neglects substitu- 
tional disorder in the Mn positions which can have a substantial impact on 
the ferromagnetism. That this is true is evident from the fact that magnetic 
properties of presently available samples are often sensitive to the conditions 
of their fabrication, and reproducibility is achieved only if the growth pa- 
rameters are carefully controlled Moreover, recent studies of post-growth 
annealed (Ga,Mn)As samples have revealed that the magnetic as well as 
the structural | |53[ | properties can depend crucially on the type of defects and 
disorder present in the system. 
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The spin- wave theory presented in the previous section describes coUective 
excitations of the ion spin system (modeled as a continuum) in terms of 
Gaussian fluctuations around the ferromagnetic ground state in the many- 
body path integral. This non-interacting spin-wave theory is exact at low 
temperatures, where deviations from the ordered ground state are small, 
but is less reliable when the temperature is raised toward the ferromagnetic 
transition temperature. 

In this section we complement the theoretical approaches described above 
by methods which (i) take disorder effects due to the randomly chosen ion 
positions into account and (ii) are able to address directly the region of larger 
deviations of the spin configuration from the ferromagnetically ordered state. 
Specifically, we present results of Monte Carlo simulations |^,^ which treat 
the minimal model without approximation. Finally we report on a rigorous 
stability analysis of the perfectly ferromagnetically ordered collinear state of 
Mn ions in the presence of disorder|^^. We predict that noncollinear fer- 
romagnetism is common in (III,Mn)V semiconductors] the robust collinear 
ferromagnetic states that have the highest ferromagnetic transition temper- 
atures occur only when the carriers are relatively weakly localized around 
individual Mn ions. 



6.1 Model considerations 

The considerations in this section are based on the kinetic-exchange model 
discussed in the previous sections. To account for the finite extent of the 
Mn ions |^ in the exchange term we replace Jpds{Ri) ■ Sj hy J d?rJ{r — 
Ri)s{r) ■ Si with the finite-range exchange parameter 

^(^) - TT^re-^^/^^'^") . (37) 
(27ra^)2 

Both the strength Jpd and range oq of this interaction are phenomenolog- 
ical parameters to be fixed by comparison with experiment or, ideally, to 
be extracted from first principles electronic structure calculations. Note that 



the exchange-coupling range parameter in Eq. (37) is required in our 
calculations once the discreteness of the Mn ions is acknowledged; exchange- 
coupling shifts of quasiparticle energies would diverge otherwise. Given this 
finite range, the only approximation we make below in treating the mini- 
mal model is that we treat the Mn spins classically. Because of the relatively 
large value of the Mn ion spins, this approximation should have minimal con- 
sequence except for the leading low temperature magnetization suppression. 
We are interested in thermal expectation values of the form 

/ = ^^ rf^^ d^?sin^Tr{/(z9,^)e-'^^} , (38) 

where (3 is the inverse temperature, Z the partition function, and i?, Lp are 
shorthand notations for the whole set of classical spin coordinates. The quan- 
tity /(i?, ip) is a function of the ion spin angles and an operator with respect 
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to the quantum mechanical carrier degrees of freedom over which the trace is 
performed. In practice we replace the fermion trace by a ground state expec- 
tation value, since the temperatures of interest will always be much smaller 
than the Fermi energy. For typical carrier densities p of order 0.1 nm~^, the 
Fermi temperature for the carriers is typically larger than lOOOK, compared 
to ferromagnetic critical temperatures lOOK. Thermal effects in the carrier 
system are therefore negligible. Thus, 

/ = 1 rdi?sin7?(0|/(i?,^)|0)e-''<*"^l°>, (39) 

-2 Jo Jo 

where |0) denotes the groundstate of non-interacting fermions with the ap- 
propriate band Hamiltonian and a Zeeman-coupling term h whose effective 
magnetic field Hcs is due to exchange interactions with the localized spins, 

h = gfiB J d^r s{r) ■ H^sir) 

H,s{r) =Y.J(''- Ri)Sni/{9fiB) (40) 

where = (sinfl/ cos 0/, sinfl/ sin (/>/, cos 0/) is the direction of the classical 
spin at Ri. In the following we denote thermal expectation values of quanti- 
ties defined in terms of classical spin orientation variables by (•) and quantum 
mechanical expectation values within the carriers ground state by (0| • |0). 



6.2 Remarks on the Monte Carlo method 

A standard way to evaluate expectation values of the form Eq. (|3^) is to 
use classical Monte Carlo algorithms which perform a random walk in phase 
space of the classical variables {"ffjip). The probabilities governing this Monte 
Carlo dynamics are specified by the dependence of many-fermion energy on 
the localized-spin configuration. The many-fermion ground state is a Slater 
determinant whose single-particle orbitals are the lowest energy eigenstates 
of a single-band or multi-band Hamiltonian. For the case of a parabolic band, 
the matrix elements of the corresponding one-particle Hamiltonian in a plane- 
wave basis read 

{k'<j'\H\ka) = ^Sk'kS.,^ + E Jk-k'e'^''-'''^'''Oi ■ t„,„ , (41) 

where k and a denote wavevector and spin indices, respectively, Jfc is the 
Fourier transform of J(r), and L the edge length of the simulation cube. 
Periodic boundary conditions restrict the admissible values of wavevector 



components to integer multiples of 2-k/L. In Eq. (41), r is the vector of Pauli 
spin matrices. 

Since the many-particle ground state of the carrier system has to be com- 
puted at each Monte Carlo step, the computational effort required for the 
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present calculations is much larger than in simple classical spin models. In 
the usual Metropolis algorithm, a single spin orientation is altered at each 
step. If this algorithm were employed here, the time required to diagonalize 
the single-particle Hamiltonian each time would severely limit the efficiency 
of the algorithm. We therefore use the Hybrid Monte Carlo algorithm, which 
was introduced in the mid 1980's in the context of lattice field theories |5^. In 
this method all classical variables are altered in one Monte Carlo step. This 
drastically reduces the number of matrix diagonalizations required to explore 
statistically important magnetic configurations. The Hybrid algorithm is a 
powerful method for Monte Carlo simulations in systems containing coupled 
classical and quantum mechanical degrees of freedom. 

A concrete Monte Carlo simulation necessarily works in a system of fi- 
nite size. In the present case, calculation of the fermion ground state also 
requires truncation of the plane-wave expansion we use for the independent 
particle wavefunctions that diagonalize the Hamiltonian We are able 

to obtain convergence with respect to this truncation only when the range 
of the exchange interaction ao is not too short, taking into account a finite 
number of plane- wave states entering the single-particle carrier Hamiltonian 
(^l]) . The importance of these finite-size effects in real and reciprocal space 
(using periodic boundary conditions), and their interplay with the regular- 
ization parameter oq, are discussed in detail in Ref. [ ^gf . In the following 
we shall leave aside such technical aspects and concentrate on the physical 
results. 

6.3 Numerical Monte Carlo Results 

In this subsection we present numerical Monte Carlo results. We concentrate 
on the spin polarizations of the Mn ions and the carriers as a function of 
temperature and address the ferromagnetic transition. We start with the 
two-band model. 

Two-band model. Fig. |lj shows typical magnetization data as a function 
of temperature. These results were obtained for a Mn ion density of Nmu — 
l.Onm"^, a carrier density p = O.lnm"'^ in a cubic simulation volume of 
V = 540nm'^, i.e., the system contains 540 Mn ions and 54 carriers. The 
effective band mass is half the bare electron mass, and the chosen exchange 
parameter is Jpd = 0.15eVnm'^. Here and in the following Monte Carlo data 
the range parameter for the exchange coupling is oq = O.lnm. 

The main panel shows the average polarization of the Mn spins. 



(42) 
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(43) 



i. e. the thermally averaged modulus of the total ion spin, along with the 
carrier magnetization, 

1 

m = — , 
pV 

which is the ensemble average of the modulus of the total ground-state car- 
rier spin. Both quantities are divided by the number of particles and are close 
to their maximum values at low temperatures. At higher temperatures they 
show the expected transition to a paramagnetic phase. The critical tempera- 
ture of this ferromagnetic transition is most readily estimated from numerical 
results for the magnetization fluctuations: 
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These two fluctuations per particle are plotted in the insets of Fig. 



3.0 



2.5 



2.0 - 



i5 1 5 - 

Mn ions 

1 1.0 - 



0.5 



0.0 




50 100 150 200 
Temperature [Kelvins] 



250 



300 



Fig. 14. Magnetization curves for Mn ions and carriers. The upper and lower inset 
show the magnetic fluctuations for the Mn ions and the carriers, respectively. Both 
differ by a factor of approximately 25 reflecting the square of the ratio of spin 
lengths. The density of Mn ions is A^Mn ~ l.Onm"^, the carrier density is p = 



O.lnm in a cubic volume of 1^ = 540nm 
mass with an exchange parameter of Jpd = 



The band mass is half the bare electron 
: 0.15eVnm^ 



They both show a pronounced peak at a temperature T ~ lOOK, defining 
the finite-system transition temperature for these model parameter values. 
In fact, in a region around this transition the two data sets differ by a factor 
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of approximately 25, which is the square of the ratio of the two spin lengths 
entering the expressions ( ^ ) and (^), respectively. This observation shows 
explicitly that the correlation length is the same for Mn ions and the carrier 
system near the transition, that is both approach the finite system size of the 
simulation. 

Our Monte Carlo approach clearly reproduces the expected ferromagnetic 
transition. Ferromagnetism still occurs when randomness in the Mn positions 
is accounted for. The transition temperature Tc can be determined unambigu- 
ously and consistently from the positions of very pronounced peaks in total 
magnetization fluctuations of both the Mn ions and the carriers. 

Results for Tc- We now turn to the transition temperature Tc for the 
two-band model. Within mean- field theory this quantity is given by Eq. (^. 
Our objective here is not to make a quantitative prediction of the critical 
temperature for particular ferromagnetic semiconductor systems. By doing a 
numerically exact calculation for a model that captures much of the physics, 
however, we hope to shed light on the range of validity and the sense and 
magnitude of likely corrections to mean-ficld-theory Tc estimates. 

The mean-field expression for Tc can be obtained by averaging the ion-spin 
and carrier polarizations over space. The effective field which each Mn spin 
experiences due to a finite carrier polarization is then constant in space and 
the carrier bands are in turn rigidly spin split by Z\ = JpdNunS- (The limit 
in which mean-field theory is exact can be achieved in our model by letting 
ao ^ oo in Eq. (|37|).) Mean- field theory, which neglects spatial fluctuations 
and correlations between carriers and Mn spins, predicts that Tc is quadratic 
in the exchange parameter Jpd and linear in the effective band mass m*, 
when correlations in the itinerant system are neglected. 

Here we determine the critical temperature Tc with the help of our Monte 
Carlo scheme. In Fig. |l^ we show results for Mn densities A'mii — l.Onm^'^ 
and a mean-field band splitting A ~ JpdNunS — 0.5eV. The left panel 
shows the dependence of the critical temperature on the carrier effective 
mass. This dependence is very important for the search for diluted magnetic 
semiconductor systems with Tc's larger than room temperature. In the mean- 
field approximation, Tc grows linear with increasing mass. The Monte Carlo 
results clearly deviate from this prediction suggesting a saturation of Tc at 
carrier masses close to the bare electron mass. For even higher masses, we 
expect the electrons to behave more classically and localize around individual 
Mn spins, suppressing long-range order further. In this limit the electronic 
energy cost of changes in the relative orientations of nearby Mn spins will get 
smaller causing Tc to decline, and eventually ferromagnetism will disappear. 
This is consistent with the observation that, within the continuum model, the 
spin stiffness declines as 1/ m* for large band masses [ ^6|j4^ , ^ . As mentioned 
earlier, the spin stiffness tends to be remain substantially larger when coupled 
light and heavy holes are retained in the calculation. Still, this calculation 



Ferromagnetism in (III,Mn)V Semiconductors 



37 




500 







0.0 



0.5 1.0 
Effective Mass m*/m 







0.1 0.2 0.3 
Carrier Density p [nm^^] 



0.4 



Fig. 15. The critical temperature Tc as a function of the carrier mass (left panel) 
and the carrier density (right panel). The exchange parameter is in both cases 
^pd = 0.2eVnm^ leading to a zero temperature mean-field spin splitting of Z\ = 
J-pdNyinS — 0.5eV. The results of the Monte Carlo runs are compared with the 
mean-field predictions. 

demonstrates that mean-field-theory must be regarded with some caution 
and its validity must be checked in each new circumstance. 

To discuss the critical temperature as a function of the exchange coupling 
parameter Jpd, we observe that the Hamiltonian of itinerant carriers satisfies 
the scaling relation 



with q > 0. Therefore the saturation of Tc as a function of the effective mass 
at fixed Jpd corresponds to a linear dependence of Tc on Jpd at fixed m*. 
This contrasts with the mean-field prediction oc J^^. 

In the right panel of Fig. |l5| we show Tc as a function of the carrier 
density. Here the Monte Carlo approach also clearly yields a lower critical 
temperatures than mean-field theory. For still higher carrier densities the 
typical distance between nearby Mn ions will become larger than the band 
electron Fermi wavelength, causing the sign of the typical exchange coupling 
to oscillate in an RKKY fashion. As we shall see in the next subsection, the 
resulting frustration makes the ferromagnetic state unstable, possibly leading 
to a regime of spin-glass order. Disordered states can also occur when the 
exchange coupling becomes strong. 




(46) 



Six-band model. We now turn to the six-band model in which the ki- 
netic energy part is given by the Kohn-Luttinger model. Fig. |l^ shows typ- 
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ical magnetization data for the (Ga,Mn)As system, assuming exchange cou- 
phng Jpd = 0.15eVnm^, carrier density is p = O.lnm^"^, and Mn ion density 
■^'ivin = l.Onm"^ in a volume oi V = 280nm''. As in the case of parabohc 
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Fig. 16. Magnetization curves for Mn ions and carriers in the six-band model for 
n exchange coupling of Jpd = O.lSeVnm^. The carrier density is p = O.lnm"^ with 
an Mn ion density of A^Mn = l.Onm^'' in a volume oi V — 280nm^. As in the case 
of parabolic bands, the ferromagnetic transition in clearly and consistently signaled 
by pronounced peaks in the magnetic fluctuations shown in the insets. 

bands, a ferromagnetic transition is clearly signaled by pronounced peaks in 
the magnetic fluctuations of both Mn ions and carriers. We find that, in con- 
trast to the parabolic two-band model, the carrier magnetization is already 
reduced at temperatures well below Tc- In fact, because of strong spin-orbit 
coupling in the valence band, full polarization of the carrier spins never oc- 
curs. Another difference compared to the parabolic-band model concerns the 
shape of the magnetic fluctuations for the Mn ions and the carriers as a 
function of temperature. Although both curves indicate the same value for 
Tc for these parameter values , their shape in the vicinity of Tc is slightly 
different and the ratio of their fluctuations is smaller than 25 (the square 
of the ratio of spin lengths involved) . These differences arise because of the 
more complicated band structure and the spin-orbit coupling present in the 
Kohn-Luttinger Hamiltonian. In the right panel of Fig. ^ we plot the transi- 
tion temperature as a function of the exchange coupling Jpd for two different 
system sizes. Both data sets agree within error bars and show a linear depen- 
dence of Tc on Jpd- This finding is the same as for the two-band model and 
contrasts with mean-field theory which predicts Tc oc J^^. The left panel of 
Fig. |l^ shows the transition temperature as a function of Jpd for the parabolic 
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Fig. 17. The right panel shows critical temperature Tc as a function of the exchange 
parameter Jpd for the same particle densities as in Fig. ^ and two different system 
sizes. Both data sets agree within error bars and show a linear dependence of Tc 
on Jpd. In the left panel, the corresponding data for the parabolic two-band model 
with an effective mass of half the bare electron mass is plotted. The latter system 
is a reasonable approximation to the six-band case. 



model for an effective mass m* = 0.5me. This value is close to the heavy-hole 
mass in the Kohn-Luttinger model for parameters appropriate for GaAs. The 
data in the left panel can be obtained from the left panel of Fig. ^ via the 
scaling relation (^6|). Comparing the two panels of Fig. |l^ demonstrates that, 
in the range of carrier densities studied here, a single parabolic band with 
an effective mass close to that of the heavy-hole Kohn-Luttinger-model band 
provides a reasonably good approximation to the behavior of the six-band 
system. We expect larger differences to occur in the strong coupling regime, 
which presents technical difficulties to the Monte Carlo calculations. 



6.4 Disorder effects and noncollinear ferromagnetism 

The magnetization data of Fig. ^ were obtained by averaging over five dif- 
ferent realizations of the Mn ions. In fact over the range of parameters we 
have explored, except for the larger values of p/Nyin, results for different dis- 
order realizations differ only very weakly from each other. This is illustrated 
in F ig. [l8| , where the magnetization curves underlying the averaged results of 
Fig. |l4| are plotted. Those five datasets are hardly distinguishable from each 
other. Positional disorder in the localized magnetic moments of Mn ions can 
however affect the nature of the ground state itself in some circumstances. 
As seen from Fig. the ion spin magnetization seems to saturate at low 
temperatures at a value slightly smaller than its maximum, given by the 
Mn spin length of 5/2. Moreover, as shown in Fig. ^ this behavior occurs 
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Fig. 18. Magnetization curves for five diflterent realizations of Mn positions under- 
lying tiie averaged data of Fig. 

quite consistently for different disorder realizations. As we see below, the 
effect is due to the randomness in Mn positions which, combined with the 
RKKY-like oscillations in the effective coupling between Mn spins, makes the 
perfectly ferromagnctically ordered collinear state state unstable and leads 
to noncollinear Jerromagnetism. 

We now outline a theory of magnetic fluctuations around a given state 
within a path integral formalism similar to the spin- wave approach presented 
in Section ^.2 for fluctuations around the collinear ferromagnetic state. In 
contrast to the spin-wave calculations, here we do not approximate the Mn 
magnetic moments by a continuum but retain them, as in the Monte Carlo 
approach, as individual spins of length S = 5/2 placed at arbitrary locations. 
Let us first consider magnetic fluctuations around the perfectly ferromagnct- 
ically ordered collinear state having all Mn spins oriented in parallel. 

Repeating the steps described in Section ||, but keeping the Mn spins 
as individual objects, one obtains the following expression for the fluctuation 
part of the effective action in up to second order in the bosonic spin variables: 

ScB = ^/ (i^m ) Dj} {v,n )zj{v„i) (47) 

where Vm — 27rm//3 is a Matsubara frequency. zi{i'm) stands for the bosonic 
Holstein-Primakoff field of Mn spin / that describes deviations from a fully 
aligned state. The fluctuation matrix Djj{i'm) reads 



Magnetisation curves for 




(48) 
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with 




(49) 



(50) 



Here (s(r)) is the expectation value of the carrier spin density, / the Fermi 
function, and 



with Tpaaif) being the spin component a of the carrier wave function with 
label a and energy = ?7q + /i, where fi is the chemical potential for the 
carriers All quantities referring to the carrier system are to be evaluated for 
the collinear orientation of Mn spins. 

The zero-frequency (to — 0) contribution to the effective action (^^ de- 
scribes the energy of static fluctuations around the collinear state. For this 
state to be stable, the matrix DYj{0) must have non-negative eigenvalues 
only, while the occurrence of negative eigenvalues of this matrix indicates 
that the perfectly collinear state is not the ground state. This interpretation 
of the zero-frequency fluctuation term is confirmed by the observation that 
this contribution is also obtained (at zero temperature) by a standard per- 
turbation theory for the carrier ground state energy with respect to small 
deviations of the ion spins from collinear orientation. The formalism given 
above embeds this finding in a more general theory of dynamic fluctuations at 
finite temperature. However, we shall concentrate in the following on static 
ground state properties, i.e. on T = 0, where the Fermi functions become 
step functions. 

We note that for any arrangement of the Mn positions i?/, the matrix 
Djj(0) contains a zero eigenvalue corresponding to a uniform rotation of 
all spins. The eigenvalues of Djj{0) are proportional to magnetic excitation 
energies. In this sense the eigenvalue distribution of Djj (0) can be interpreted 
as a density of states for magnetic excitations. 

We have evaluated the spectrum of Djj{0) in systems given by a simula- 
tion cube with periodic boundary conditions averaging over different realiza- 
tions of the Mn positions. The single-particle wavefunctions ipaa-i'i') are com- 
puted in a plane-wave basis taking into account wavevectors q with length up 
to an appropriate cutoff qc. The same truncated plane- wave basis is used to 
compute the quantities ( |5l| ) entering (pO|). Note that, for fluctuations around 
the collinear ferromagnetic state, DJj (ioj) is always real and symmetric for 
real to since all carrier wavefunctions have, for a given spin projection a, 
the same coordinate-independent phase. This follows from the fact that the 
single-particle Hamiltonian describes for each spin projection just the prob- 
lem of a spinless particle in a potential landscape provided by the Mn ions. 




(51) 
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Since Djj{iuj) is real and symmetric, the components of each of its eigen- 
vectors all have the same phase (and can be chosen to be real). Physically 
this corresponds to the invariance of the system under rotations around the 
magnetization axis of the coUinear state. 

The two upper panels of Fig. p^show results for typical system parameters 
for two different values of Qc- Comparison of the panels shows that the effects 
of the wavevector cutoff on the low-lying excitations have already saturated 
for the smaller Qc- Almost all eigenvalues of Djj{0) lie at positive energies, but 
there is often a small fraction of negative eigenvalues, indicating an instability 
of the coUinear state. 

In the calculations discussed so far, the Mn positions were chosen com- 
pletely at random with a uniform distribution, while in a real (III,Mn)V 
semiconductor the Mn ions are supposed to be located on the cation sites 
forming an fee lattice. In the bottom panel of Fig. |l^ we show data for the 
same system parameters as in the top panel but with the Mn positions chosen 
from an appropriate fee lattice such that about 5% of all sites are occupied. 
Both plots are practically identical indicating that our observations do not 
depend on this detail of the modeling. 

The shape of the eigenvalue distribution of the fluctuation matrix Djj (0) 
is, in the model we have studied, sensitive to the Mn density Nmu, the car- 
rier density p, and the Hamiltonian parameters m*, Jpd-, and ag. Situations 
in which the coUinear ferromagnetic state is, for certain disorder realization, 
stable can be approached most simply, for technical reasons, by letting ao be 
larger. However for the value of p/Nmh illustrated in Fig.|l^, corresponding 
to a carrier density somewhat lower than measured in the highest Tc sam- 
ples, negative eigenvalues occur for nearly any Mn ion distribution. We note 
that negative eigenvalues increase in number as the wavevector cutoff is in- 
creased toward its converged value. These results suggest that noncoUinear 
ferromagnetic states are common, that they are sensitive to the distribution 
of Mn ions and other defects - especially those that trap carriers, and that 
the coUinear ferromagnetic state tends to become unstable as mean-field band 
eigenfunctions become more strongly localized around Mn ion site 

To further analyze the nature of this instability we consider the partici- 
pation ratios for these excitations which we define by 



(52) 



where a j is the /-th component of the j-th normalized eigenvector of Djj (0) . 
The ratio pj is an estimate for the fraction of Mn sites that have important 
involvement in the j-th spin wave. For example the zero-mode, uniform rota- 
tion of all spins, haspj = 1. Fig.pO| shows the disorder-averaged participation 
ratio as a function of spin- wave energy for the same situation as in Fig.p^. 
The property that negative energy excitations have large participation ratios 
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Fig. 19. The disorder-averaged eigenvalue density of the matrix Djj{0) describing 
magnetic fluctuations around the coUinear state. The data are obtained for a sim- 
ulation cube of volume V — L"^ = 400nm'^ with a Mn density of A'Mn ~ l.Onm"^ 
and a density of p = 0.15nm~^ of carriers having a band mass of half the bare elec- 
tron mass. The strength of the exchange interaction between ions and carriers is 
Jpd = 0.05eVnm~"^ with a spatial range of ao — 0.40nm. The two upper panels show 
data for different wavevector cutoff qc with the Mn positions chosen completely at 
random. The lowest panel contains data for the same situation as the top one but 
with the Mn positions chosen from an fee lattice. The peaks at zero energy are due 
to the uniform rotation mode which strictly occurs in any disorder realization. 



shows that the instabilities of the coUinear state involve correlated reorienta- 
tions of many spins, rather than lone loosely coupled moments. 

Fig. ^ shows the disorder-averaged participation ratio for the same situ- 
ation as in the top panel of Fig. |l9|. The negative-energy modes have clearly 
higher participation ratio than the eigenvectors at positive energy. This shows 
that the instability of coUinear state is due to long-ranged fluctuations involv- 
ing a large fraction of the spins present in the system. Qualitatively the same 
observations are made for other values of system parameters. 

The effect of a weak external magnetic field on the spin-wave excitation 
spectrum of the noncoUinear state is particularly simple. The field couples to 
the local moment through its Lande g-factor, adding ifiBHext to the energy 
of a spin- wave excitation for S — 5/2, and to the spin and orbital degrees 
of freedom of the band electrons. The orbital coupling leads to Landau lev- 
els that do not play an important role at weak fields in these highly disor- 
dered samples. Zeeman coupling is also unimportant, yielding a contribution 
to the spin-splitting that is negligible compared to the mean-field splitting 
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Fig. 20. The disorder- averaged participation ratio for tlie same situation as in the 
top panel of Fig. ^ The data is averaged over the sample intervals of the histogram. 
The value at zero energy is enhanced due to the contribution of the uniform rotation 
mode in that sample interval. 

^ = JpdNMnS ~ O.leV. It follows that the most negative eigenvalue of 
Djj is the value of g^idBHext necessary to force full spin alignment of a non- 
collinear state. The experiments of Potashnik et al. demonstrate that the 
maximum value of Af (T = 0) is achieved over a certain range of annealing his- 
tories. We associate this maximum value with the fully aligned collinear Mn 
configuration state; indeed the maximum moment per Mn is consistent with 
full alignment partially compensated by band electrons. For other annealing 
histories, M{T = 0) is reduced, corresponding to noncollinear order of the Mn 
spins. Calculations like those described above show that these spins gradually 
align as an external field is added to the Hamiltonian. We expect full align- 
ment to be indicated experimentally by a kink in the M (T — 0, Hext) curve. 
At this point, we predict that the system will still have gapless excitations and 
power-law temperature dependence of the magnetization, in sharp contrast 
to the gapful excitations and exponentially suppressed temperature depen- 
dence that holds for conventional ferromagnets in an external magnetic field. 
Spin-resonance experiments will nevertheless see a gaped spin- wave spectrum 
in this regime, since they couple only to the zero mode. 

Finally we comment on the degree of spin alignment in the noncollinear 
state. We generalize our formalism by expanding around and defining Holstein- 
Primakoff bosons with respect to a spin coherent state configuration with 
orientations i7/ = (sini?/ cos 95/, sin??/ saiLpi, cosi?/). In this general case the 
effective action has a static contribution (at zero Matsubara frequency only) 
linear in the Holstein-Primakoff variables: 




(53) 



with gi = g] + igj 



j, and 
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{{{sir)) ■ + ((s(r)) •e,)e,) X (54) 

g] = V2Se, ■ (^e^, x J d\j{r ~ Ri){s{r))^ (55) 

where e^^ — (cos ip i, sin (pi, 0) and — (0, 0, 1). 

The components gi represent the gradient of the energy with respect 
to distortions parametrized by the zj. As above in the case of fluctuations 
around the coUinear state, this zero-frequency contribution can also be ob- 
tained via perturbation theory around the given state of Mn spin orienta- 
tions. The contribution to 5fluc quadratic in the Holstein-Primakoff variables 
involves again a proper Matsubara sum and can be obtained straightfor- 
wardly. However, its concrete form for the general case is more complicated, 
and shall not be analyzed here. A given orientation of Mn spins is stationary 
with respect to fluctuations if all complex coefhcients gj vanish. This is the 
case if and only if J (PrJ{r — i2/)(s(r)) is coUinear with the direction Hj 
of the local ion spin. Therefore, the coUinear ferromagnetic state is always a 
stationary (but not necessarily stable) spin state. 

We have employed the energy gradient expression ( |53|) in a numerical 
steepest descent procedure to search for true energy minima. Our results are 
as follows: In cases where the energy minimum found by this method is close 
to the coUinear state (with a magnetization of about 90% of the maximum 
value or more), this minimum appears to be unique for each disorder re- 
alization. We can therefore be confident that we have located the absolute 
ground state of the system. In situations where the magnetization is reduced 
more substantially, however, (by ~ 20% or more) we converge to different 
energy minima from different starting points. In these cases the model has 
substantial spin-glass character with a complex energy landscape. This situa- 
tion occurs typically for larger ratios p/A^A/n- For the system shown in Fig. 
for instance, magnetization values at local energy minima are typically 30 to 
40% of the coUinear state value. 



7 Concluding remarks 

III-V compound semiconductors with Mn substituted for a small fraction of 
the cations are ferromagnets. In this chapter we have primarily discussed the 
picture of the physical properties of these materials that follows from a model 
in which each Mn produces a S=5/2 local moment and acts as an acceptor. 
The valence band holes and the exchange interaction that couples them to 
the local moments are both treated phenomenologically. In the simplest ap- 
proximation the hole bands are taken to be those of the host semiconduc- 
tor, so that this model has a single unknown parameter, the exchange cou- 
pling strength, whose value can be determined experimentally. Many physical 
properties of this model are relatively simply calculated when a virtual crys- 
tal approximation is adopted. When an envelope function approach is used 
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to describe the semiconductor bands this is equivalent to replacing the Mn 

ions by a continuum, completely eliminating disorder. Wc have shown that 
the mean-field theory of this model predicts critical temperatures in good 
agreement with experiment, that it can describe qualitative features of the 
magnetic anisotropy energy including the effect of lattice-matching strains 
between the magnetic epilayer and its substrate, and that it implies values of 
the anamolous Hall effect often used to characterize these ferromagnets that 
are in good agreement with experiment. We have also described a theory of 
the elementary collective spin-wave excitations of the virtual crystal model. 
By comparing the energy of these elementary excitations with those of the 
mean-field theory, we are able to judge the reliability of mean-field-theory 
critical temperature estimates. In this way we find that we reach the wrong 
conclusion about mean-field theory when we make the apparently innocent 
choice of using a single parabolic band whose mass is equal to the heavy hole 
mass of the host semiconductor. For such a model, we find that with typical 
parameters mean-field theory would give a rather poor estimate of the crit- 
ical temperature. However when we use a realistic band model with heavy 
and light holes that are coupled the typical energy of collective magnetic 
excitations increases substantially. The fact that heavy hole states are given 
approximately by J=3 /2 spin coherent states with a k dependent orientation, 
plays an essential role in producing the enhanced spin stiffness. We conclude 
that a realistic treatment of the valence band, including its characteristic and 
strong spin-orbit mixing, is absolutely necessary to achieve even a qualitative 
understanding of critical temperature trends in these ferromagnets. 

Some properties of (ni,Mn)V ferromagnets are misrepresented by the 
virtual crystal approximation. We have examined the effect of Mn site disor- 
der on these ferromagnets both by performing Monte Carlo calculations and 
by evaluating the collective excitations of disordered systems. The Monte 
Carlo calculations evaluate the properties of the model exactly for finite sys- 
tems, except for treating the Mn spin orientations as classical which will have 
a small effect on thermal magnetization suppression at low temperatures. 
These calculations confirm conclusions reached by using the virtual crystal 
approximation for the most part, but also highlight some of its limitations. 
In particular, while the coUinear magnetic state is always stable in the vir- 
tual crystal approximation, it can be unstable when disorder is accounted for, 
leading to non-coUincar but still ferromagnetic states. This property provides 
an explanation to the dependence of a given sample's transition temperature 
and of the magnetization on details of growth and post-growth annealing 
procedures. 

Research on Mn-doped semiconductors is a very active area of physics, 
both theoretically and experimentally and there are many topics which we 
have not been able to even touch upon. For example, studies of nano structured 
semiconductor systems with Mn-doped components represent an important 
component of developments that are aimed towards future spintronic de- 
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vices. Nanostructures such as magnetic quantum wells are expected to be 
unusual ferromagnets because of the possibility of using confinement effects 
and doping profiles to manipulate their magnetic properties. The transition 
temperature of these systems can be tuned by external electric field p7| , as 



demonstrated in recent studies of (In,Mn)As field effect transistors [58 . It is 



also predicted |59| that hysteresis properties of magnetic quantum wells will 
be extremely sensitive to external bias voltages. As a result, the magneti- 
zation orientation in quantum wells can be manipulated electrically without 
changing the magnetic field. The possibilities for nano-engineering of ma- 
terial properties have already been made apparent by the relatively simple 
(Ga,Mn)As digital ferromagnetic heterostructures l^^l- These systems are 
grown by incorporating submonolayer planes of MnAs into a GaAs host. 
They show ferromagnetic transition temperatures up to 50 K ||6l| , |62| , the 
anomalous Hall effect and remarkably square hysteresis loops with higher 
coercivities ||6l|| than their random alloy counterparts. A confident modeling 
of these alternative ferromagnetic semiconductors with controlled Mn distri- 
bution would significantly contribute to our understanding of magnetic and 
transport properties of (III,Mn)V ferromagnets. 

In this Chapter we have highlighted one particular approach to modeling 
these materials. We anticipate that other approaches will also bring useful 
insights. In particular dynamic mean-field-theory calculations | |30t are able 
to capture some of the Mn alloy disorder physics in a relatively simple way, 
although they do not capture spatial correlations and would not capture 
non-coUinear states. Similarly first principles calculations |]l8j are likely to 
prove extremely valuable in the future, once LDA-I-U or dynamical-mean- 
field theory corrections have been applied to treat the Mn d-electron local 
moments more accurately. 
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